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Entropy generation in quantum sytems is tied to the existence of a nonclassical 
environment (heat bath or other) with which the system interacts. The continuous 
'measuring' of the open system by its environment induces decoherence of its wave 
function and entropy increase. Examples of nonrelativistic quantum Brownian 
motion and of interacting scalar fields illustrate these general concepts. It is shown 
that the Hartree-Fock approximation around the bare (h = 0) classical limit can 
lead to spurious semiquantum chaos, which may affect the determination of entropy 
production and thermalization also in other cases. 

1 Introduction 

The related topics of environment induced quantum decoherence and the tran- 
sition from quantum to classical mechanics have recently been investigated in 
widely varying contexts, ranging from problems of interpretation of quantum 
mechanics, of the measurement process in particular, to questions of how clas- 
sical physical laws and the observed classical features, of cosmology emerge 
in the underlying quantum universe, see e.g. Refs.EJ-0. The central idea of 
the quantum decoherence approach is that the transition quantum — > classi- 
cal can be understood as a dynamical effect within quantum mechanics itself, 
especially keeping h nonzero as it is. 

In this lecture I will explain some of these issues, in particular how entropy 
generation is intimately related to an open quantum system. Its identification, 
its properties, and the calculation of the time dependent entropy production 
present unresolved problems for strongly interacting systems, such as the high 
energy density matter formed during relativistic hadronic or nuclear collisions. 
Based on more technical details given in Refs.Q, I will outline how the quantum 
decoherence approach can be applied here. 

This is obviously related to the persistent issue of thermalization of strongly 
interacting matter, the assumption of which is one of the conceptual corner^ 
stones guiding the ongoing search for the quark-gluon plasma (QGP), see Rcf.Q 
and earlier "Quark Matter" proceedings. As reflected during this conference, 
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many theoretical attempts to pinpoint the most relevant QGP properties de- 
part from the assumption of a thermalized high entropy density state of matter. 
However, at present we are still quite far from understanding how an initially 
'pure' quantum state with zero entropy, e.g. two colliding nuclear wave packets, 
could evolve into the 'mixed' QGP state. 

In order to fully appreciate the problem, we need a precise definition of 
the entropy, which preferably should allow us to extrapolate to the usual ther- 
modynamical limit or Boltzmann's entropy of statistical mechanics. Let us 
recall that the First Law of Thermodynamics relates infinitesimal changes of 
the internal energy U, the volume V, and the entropy S, given the pressure P 
and the temperature T of the system: 



Thus, the change of the internal energy is related to the work done by the 
system against the external pressure and to the heat 5Q transfer or the en- 
tropy change, respectively. The Second Law of Thermodynamics rules that 
the entropy in a closed system cannot decrease: 



Furthermore, based on a statistical definition of the (TV-particle) entropy, 



employing a suitable phase-space distribution function, one is led to the in- 
terpretation due to Boltzmann that the entropy measures the observer's lack 
of information about the system. It determines the number of possible real- 
izations (microstates) of the system which conform with a given macrostate 
specified, for example, in terms of the above thermodynamical variables. 

Evidently, in relativistic heavy-ion physics much use is made of such re- 
lated thermodynamical notions as equilibration, thermal pj_-spectra, equation 
of state, phase transition, etc. 

It is our aim here to get a quantum mechanical handle on the entropy. - 
As is well known, the Wigner function is what comes closest Jjo a quantum 
mechanical distribution function replacing f(x,p;t)a; see RefJ3 for a review 
of earlier work on quark-gluon transport theory, i.e. the dynamics of QCD 
Wigner functions. For our present purposes it is more convenient to work 
with density matrices, which generally are related to Wigner functions by ap- 
propriate Fourier transforms. - In terms of the density operator p character- 
izing a (closed) quantum mechanical system, we will henceforth employ the 



dU = -PdV + SQ = -PdV + TdS . 
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von Neumann definition of the entropy 



01: 



S(t) — — Tr p{t) lnp(t) = — w n lnm, 



n 



(4) 
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where the trace is over a complete set of states and w n denotes the probability 
of such a state \n) ; from now on we employ units such that h = c = ks = 1 , 
except when stated otherwise. 

As we shall see in the next section, the von Neumann entropy has several 
desirable properties, but shows some surprising features as well. Most notably 
S =JX for closed systems. The latter result has led to some confusion in the 
pastE2l, in particular it has been claimed that S, as defined here, is completely 
useless, if one wants to characterize the entropy production in high-energy 
physics. However, these objections are overcome by the decoherence approach 
and especially by realizing that in this case, like in most if not all physically 
interesting cases, the system is open indeed. 

2 The von Neumann Entropy and the Need for Open Systems 

In order to illustrate some important features of the entropy S, as defined 
in eq. (Q), we evaluate it for two limiting cases, for a quantum system in a 
pure state and in the thermal equilibrium state, respectively. An elementary 
example for a two-state system is given in Ref B(c) . 

If the system is in a pure state \^), the density operator assumes the simple 
form p = \$>)(?$>\. Correspondingly the probability of finding the system in 
this state is = 1. Therefore, we obtain immediately: 



which can be interpreted that we know everything about the system that we 
possibly can. 

Similarly one finds that if the system is in any one of N states with equal 
probability l/N, then the entropy equals IniV , i.e. the maximum value, when 
we are completely ignorant about which state the system is in. 

Let us consider one of the most important cases of an impure or mixed 
state, namely when the system is in thermal equilibrium. Then, the density 
operator is given by: 



5* = -1-lnl = , 



(5) 
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where Z is the partition function, H is the Hamiltonian, and (3 = T 1 denotes 
the inverse temperature of the system. We find: 

_ e -0E n -f3E n -0H 

5 ( T ) = -£ ^T ln ^T = ^Z-Tr—{-[iH) 

n 

= \nZ - Z- l j3dpZ = d T (T\nZ) = -drF , (7) 

which relates the von Neumann entropy to the usual partial derivative of the 
free energy, which is familiar from thermodynamics. 

The properties of the entropy which we have considered so far tie in with 
our knowledge of thermodynamics or classical statistical mechanics. However, 
a problem arises immediately, since the entropy S turns out to be a constant 
of motion. According to the Schrodinger equation the time evolution of wave 
functions is unitary, | <!>(£)) = exp(—iHt)\$>(0)) , which implies: 

S(t) = -Tr {c- l6t p(0)e lAt \n[e- l6t p(0)e lAt } = S(0) , (8) 

using the cyclic property of the trace. In a closed system the entropy stays 
constant at the initial value. 

More specifically, the unitary quantum mechanical time evolution prohibits 
a transition from a zero entropy pure initial state to a mixed nonzero entropy 
final state. However, it is well known how to overcome this impasse. Similarly 
as in classical statistical mechanics, where coarse graining in phase space allows 
to circumvent df/dt — (Liouville flow) and hence Sb = 0, cf. eq. (H), we 
have to give up the premise of unitary time evolution. 

The decoherence approach is based on a separation of 'all degrees of free- 
dom' into the observed degrees of freedom (system) and the 'rest of the uni- 
verse' (environment) EJoD. The border between system and environment gen- 
erally has to be open to the exchange of information and possibly to the ex- 
change of energy, momentum, etc. Without interaction two isolated closed 
systems result. Then, the density operator p representing all degrees of free- 
dom (d.o.f.) evolves unitarily and the corresponding entropy is a constant, as 
we have shown. However, the density operator p$ representing the system, 

ps(t) = Tr EP (t) , (9) 

which is obtained by tracing over the environment degrees of freedom, evolves 
in a more complicated non-unitary way. Consequently, the entropy of the 
system Ss(t) , which is defined as the von Neumann entropy employing ps(t) , 
will change in time, as we shall see. 
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The non-unitary evolution of the system density matrix is accompanied by 
environment induced decoherence, i.e. typically its off-diagonal elements which 
encode the quantum mechanical interference effects decay with a character- 
istic decay constant to ■ Correspondingly, the eigenvalues of ps change and 
the diagonal matrix elements become probabilities characterizing the resulting 
mixed state. If there is thermalization, 



pT(t) ^ » (io) 



i.e. the diagonal matrix elements become the usual Boltzmann weights, cf. 
eqs. @, (0). Most interestingly, the decoherence time Td is found to be several 
ten orders of magnitude smaller than the thermalization time for models of 
macroscopic bodies interacting with theix-environment (cosmic background ra- 
diation, atmosphaeric gas particles, etc. til. This explains the classical behavior 
of most of what we observe in daily life u. Cases of quantum decoherence and 
'revival' of the off-diagonal matrix elements are also known and experimentally 
observed in cavity quantum electrodynamics. 

For microscopic systems, generally, the situation is more intricate and 
depends sensitively on the interactions and the separation into system and en- 
vironment d.o.f., in particular for more complex-systems than a single particle 
or one d.o.f. interacting with an environment til. If we find td *C t*/i , then 
the possibility arises of creating a system with high entropy which is not at all 
thermalized. We would like to know precisely how these time scales compare 
in high-energy collisions of hadrons or nuclei and, of course, to what extent 
these systems thermalizc. 

Next, we present the useful theorem that the entropy generated in a system 
exactly equals the entropy generated in its environment due to the mutual 
decoherence process, 

S s {t) = -Tr p s [t) \n Ps (t) = -Tr p E {t) hxp E (t) = S E (t) , (11) 

with pe calculated from p by tracing over the system d.o.f. analogously to 
PS , eq. (Hi . This result follows from the Schmidt decomposition of the density 
matricesEl(a), which is for matrices (e.g. in Hilbert space) what the Schmidt 
orthogonalization procedure is for vectors. 

This leads us to the question, whether the hard as well as the large number 
of soft photons generated by bremsstrahlung from quarks during the initial 
hard scattering and stopping phase of a hadronic collision yield an important 
decoherence effect. This would imply that a major part of the entropy observed 
in the final hadronic states of the reaction is already generated during its 
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initiation. A simple counting argument based on the ratio of effective photon 
and QGP d.o.f., i.e. their entropy ratio in thermal equilibrium, 

N QGP 37aVT3 ' 1 ' 

seems to indicate that only about 5% (or less) of the entropy can be gener- 
ated in this way. Note, however, the strong temperature dependence of the 
absolute numbers. Even if it does not maximize the entropy, as a thermal 
distribution does, a nonequilibrium distribution generated by bremsstrahlung 
processes which effectively corresponds to a higher photon 'temperature' may 
lead to a sizeable correction of the above estimate. The discussion following 
eq. (^) explains that generally a flat distribution contains less information/more 
entropy than a steep one and that the size of the available phase space matters 
('lnJV'). A dynamical calculation generalizing results presented in the following 
Sections 3 and 4 is presently carried out. 

In the example just considered the separation between system and environ- 
ment d.o.f. is based on the fact that the strongly interacting system of quarks 
and gluons is essentially bound for, say, lOfm/c until 'freeze-out', whereas 
the nonequilibrium environment photons are generated particularly during the 
first 0.1 ... 0.5fm/c and essentially free to leave. Thus the 'rest of the universe' 
or environment is open to interpretation and has to be determined according 
to the physical circumstances. 



3 A Nonrelativistic Quark in Brownian Motion 

The purpose of this section is to demonstrate how the general concepts of 
environment induced decoherence and entropy generation in an open quantum 
system work in practice. The example can be understood as a nonrelativistic 
(heavy) quark moving and interacting with its own gluonic environment field, 
the properties of which can only be modelled at present. 

We start with a model Hamiltonian describing the bilinear translation 
invariant interaction between a nonrelativistic particle (system) and an infinite 
set of harmonic oscillators (environment): 

^ = !^+E{lt + M^-*) 2 } ' (13) 

where fl presents a high-frequency cut-off (such as Debye frequency and inverse 
classical electron radius in similar models related to the polaron and electron, 
respectively); we may think of the QCD scale parameter Aqcd ~ 200 MeV. 
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This type of model associated with the names of Feynman and Vernon or 
Caldeira and Leggett has been extensively studied in the field of quantum 
Brownian motionEj, developing further the original Feynman- Vernon influence 
functional technique. 

As it turns out, the physics described by the model of eq. ( p"3| ) is completely 
determined by the spectral density distribution of the environment: 

7 M = \ £w„<fi m n u^S(u} - uj n ) , (14) 

and the related noise and dissipation kernelsS. - For the polaron and electron 
cases one obtains I(u>) oc uj 3 for uj — » , which lead to the well known mass 
renormalization for t — > oo , whereas for an 'Ohmic' environment one chooses 
I(u>) oc w, leading to a dissipative force proportional to the velocity of the 
Brownian particle. 

Unfortunately, in the QCD case we do not know how to calculate the 
spectral density distribution. Classically QCD is a nonintegrable theory, which 
implies that such a decoupled oscillator representation of the dynamics does 
not exist E3. Therefore, it seems unlikely that standard perturbation theory, 
which is based on the ?i-(loop-)expansion around the naive classical limit, can 
provide reliable results in this context. 

Keeping this limitation in mind, I considered this model in the strong- 
coupling limit for t <C at T = , where only certain moments of I (to) enter 
the calculation and can be chosen as parameters of the model. Presently the 
spectral density is much larger in the infrared than in the polaron/electron case. 
Technical details of the calculation of the reduced density matrix propagator 
in the environment corresponding to eq. (13) can be found in Ref.El(a). Also 



longer time scales may be interesting to study. 

Here we want to discuss only the very initial influence of the environment. 
Choosing I(u>) = g£l 3 F(u> /fi)0(fi — uj) , with an undetermined shape function 
F and dimensionless coupling constant g , we have go = g f Q dx x F(x) . 
Then, the following parameter and dimensionless time variable, respectively, 
are needed: 

a=— — (w n)- i . t+ = nt(2g — f' 2 , (15) 

where wq denotes the initial (t = 0) width of a Gaussian wave packet repre- 
senting a particle moving with the velocity vq = (p)o/M initially. 

In Fig. 2 a) we observe the dissipative 'friction' effect of the environment. 
In the present approximation we obtain a universal curve showing the decel- 
eration of the particle in terms of the time dependent velocity (divided by iio) 
as a function of t+ . For g = 0, of course, v(t + ) — vq . Secondly, as shown 



7 



a) Environment Friction 




t + 



b) Wave Packet Spreading/Squeezing 




' ' ' 0.2 ' ' o'.t ' ' 0.6 ' ' 0.8 ' ' ' 1 



t + 

Figure 1: Environment effects in nonrelativistic quantum Brownian motion, a) The cen- 
ter of a Gaussian wave packet slows down, b) The width of the wave packet (a = 
3.0, 2.0, 1.0, 0.3, 0.1 , full lines, top to bottom) expands more rapidly or may be squeezed as 
compared to the usual case (a = 1.0, 0.3, upper and lower dashed curve, respectively). See 
main text for details. 
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a) Wave Paket Spread/Squeeze 
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Figure 2: (a) Similar as Fig. 2 b) , however, more extreme initial widths are considered — 
small, medium, and large wave packets (top to bottom), respectively, (b) The corresponding 
single particle entropy for a large, small, and medium size wave packet (top to bottom). See 
main text for details. 
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in Fig. 2 b), we find some surprising effects which the environment can have 
on the width of the wavepacket. There, the relative change of the width as 
a function of t + is represented: the wave packet spreads faster (a > 1) than 
a corresponding free (g = 0) quantum mechanical wave packet, the width re- 
mains constant (a = 1), or it even gets squeezed (a < 1), depending on the 
parameter a . The dashed curves show the spreading for g = 0; they do not co- 
incide because of the rescaled time axis. If the other physical model parameters 
were fixed, this result clearly shows that the particular coordinate-coordinate 
coupling between the particle and its environment, which we have chosen in 
cq. (|l3|), is sensitive to the size of the particles' wave packet: sufficiently small 
packets spread faster than usual, sufficiently large ones get squeezed. 

This is illustrated once more for some more extreme cases in Fig. 3 a), 
where the three curves (top to bottom) correspond to an initially small (wqQ <C 
1), medium size (wofl « 1), and large (wofl 3> 1) wave packet, respectively. 
For these three cases we show in Fig. 3 b) the corresponding entropy, i.e. the 
entropy evaluated with the help of the time dependent density matrix of the 
nonrelativistic particle. In this figure the curves correspond to an initially large 
(top), small (middle), and medium size (bottom) wave packet, respectively. We 
do not attribute too much importance to the sizeable entropy which can be 
reached in very short time, e.g. O.lfm/c, since the numbers depend on the 
model parameters; we recall, however, the ultrarelativistic/high-T limit with 4 
units of entropy per particle. - Actually represented here is the simpler linear 
entropy,-Su n = Tr[p — p 2 } , which provides a lower bound for the von Neumann 
entropy 0(a). The linear and von Neumann entropies coincide in the zero and 
maximal entropy limits. - It is remarkable that medium size wave packets, the 
width of which is least affected by the environment, also generate the lowest 
entropy increase. They essentially behave like a classical particle. 

Furthermore, coherent superpositions of such states decohere most quickly 
and show the strongest entropy generation u(a). These properties taken to- 
gether qualify the Gaussian wave packets as approximate pointer states of the 
present model, similarly to the simple rnndel studied in Refo. Pointer states 
show the behavior of a classical pointer EoEl: coherent superpositions of such 
states decohere very rapidly into mixed states, which are observed as various 
pointer positions with certain probabilities in a (macroscopic) experiment. 

4 Interacting Scalar Fields 

In this section we study two interacting scalar quantum fields, representing 
the system ('1') and environment ('2'), respectively. We want to describe their 
dynamics as well as the entropy generated in the system. We sketch this first 



10 



step in generalizing the Brownian motion example of Section 3 towards more 
realistic cases related to high-energy collisions, while the technical details may 
be found in Refs.E. 

Let the classical action of this model be given as: 

S = J d A x{L x + L 12 + L 2 } , (16) 

with the interaction potential L\i = — V($i,$2)j an d the Lagrangians: 

L j = $d lt $ j d»$ j -v j ($ j ) , «(<&) = + iA$ 4 . (17) 

Consider the separation into system and environment as originating from split- 
ting a single scalar field, <E> = $i + $2, according to some 'coarse graining 
prescription' to be discussed shortly. Then, the interaction potential is: 

F($!,$ 2 ) = -M 2 $i$2 -^$i9 /i $2 + TjA$i$ 2 [4$ 1 2 + 6$i$2 + 4$2 2 ] ■ (18) 

We see that the bilinear coupling studied in most quantum Brownian motion 
models appears naturally here among other terms. Presently we will keep only 
this first term for simplicity, while results for a general quartic V($i,<f>2) are 
derived in RefJj(a). - Now, various interpretations of the fields $1, $2 are 
possible: 

• Coordinate space: $1 and $2 represent the same field, however, with 
support either inside the system or in the environment, respectively; 
V($i,$2) presents a contact interaction at the surface. Models of this 
type have been studied for the geometric entanglement entropy generated 
by event horizons (black holes), see e.g. Refs.Ej. 

• Momentum space: $1 and $2 represent the long- and short-wavelength 
components of the same field. This has been ofomuch interest in cos- 
mology in the context of inflationary models; Ref.ll3 is a recent example, 
from which other references can be traced. 

Strongly interacting matter naturally leads to a separation of d.o.f. In a 
finite size QGP droplet, with radius R on the order of a few fm , quark-gluon 
modes with wavelengths much larger than R are suppressed due to confinement 
and form an environment of virtual fluctuations together with composite meson 
(baryon) fields in the nonperturbative vacuum. Its properties, as with the 
case discussed in Section 3 , cannot be calculated at present, but presumably 
play a crucial role in multiparticle production processes and the associated 
entropy generation. Finally, we mentioned before the bremsstrahlung photon 
environment which contributes to the decoherence of the charged quark system. 

In the following we assume $1 and $2 to be two physically distinct fields. 
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4-1 Variational Approach to Time Evolution in Field Theory 

We now turn to the dynamics of the coupled fields of our model defined in 
cqs. (pf)-([j"8|). The standard approach is motivated by perturbation theory, i.e. 
the usual ?i-(loop-)expansion, and based on the Schwinger-Keldysh formalism. 
In order to resum the large class of (iterated bubble) diagrams corresponding 
to the time dependent Hartree-Fock approximation (TDHF), it is much more 
efficient to employ Dirac's time dependent variational principle. 

The starting point is the functional Schrddinger equatum describing the 
full dynamics of a generic field <p in the Schrodinger picture til: 

id?B[<p;t] = H[it,<p]*[<p;t]= j 'd d z{-i^ + |(V^) 2 +V(^)}^;i] , 

(19) 

where \E'[<^;i] = (<p\i&(t)) denotes the wave functional in the (^-representation, 
which corresponds to a wave function ip(x,t) = (x\ip(t)) for a one-dimensional 
quantum mechanical system, and 7r = —id /Sip represents the canonical mo- 
mentum operator conjugate to the field ('coordinate') ip . The dynamics is 
determined by the Hamiltonian H. In this context the completeness and inner 
product relation, respectively, involve functional integrals instead of ordinary 
ones (orthogonality needs a S- functional) , for example: 



<*i(t)|¥ 2 (t)> = / £V (*iWI^^I*2(i)) = / Vip ^l[p;t}^ 2 [ip;t] , (20) 

with details to be found in Ref.i(a). 

Then, we may state the variational principle: 

<5T[¥] 



6V 



, for all * with (\&(t)|*(t)) = 1 , (21) 



where r[*] = j dt (${t)\[id t - H\\V(t)) , (22) 

i.e. requiring the effective action T defined in eq. ( p2| ) to be stationary against 
arbitrary variations of the normalized wave functional \& , which vanish at 
t — > ±oo, is equivalent to the exact functional Schrodinger equation, eq. (19) 
above. With the variational principle in hand, one can solve the time evo- 
lution problem approximately by choosing a suitably parametrized trial wave 
functional. We remark that the effective action is real by construction. 

We work with the most general Gaussian trial wave functionals. For a 
generic field (p it is defined by: 

9 a [<p;t] = N{t)exp{-[<p-<p{t)] [JG" 1 ^) - *E(t)] [<p-<p{t)}+m{t)[ V -${t)}} , 

(23) 
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where henceforth we do not explicitly write the integrations over spatial vari- 
ables. For example, 



<pG-\t)<p(t) = J & d xd d yv{x)G- x {x,y,t)(p{y,t) , 

n(t)(f = J d d x Tt(x,t)(p(x) , tr£(i) = J d d x T,(x, x, t) . (24) 
The normalization factor N (for symmetric and positive-definite G ) is: 

1= /2tytf&tot]¥ G tot] — > N(t) = (M dot Git))- 1 ^ , (25) 



cf. eq. ( p0[ ) (A/" is an infinite constant which we omit henceforth). 

The meaning of the variational parameter functions (p, n, G, and £ is dis- 
cussed in Ref.Q(a). The choice of Gaussian trial wave functionals is dictated by 
the need to evaluate functional integrals, in order to calculate the effective ac- 
tion r , eq. (^2|). The equivalence of this Hartree-Fock effective action with the 
Cornwall- Jackiw-Tomboulis generating functional for two-particle irreducible 
graphs was demonstrated in Ref.lia for energy eigenstates of the field. Variation 
of r w.r.t. to the parameter functions will give the dynamical equations, cf. 
Section 4.3 below, representing the field theory in terms of coupled equations 
for the one- and two-point Wightman functions. 

We proceed with the ansatz for the two-field wave functional: 

*ia[$i,$ 2 ;t] = iVia(*)*Gj*i;*]*G a [*2;*] (26) 
• exp{-|[#i-f 1 (*)][G 12 (t)-iE 12 (t)][*a-*2(t)]} , 



with the normalized Gaussians on the r.h.s. as defined in eq. (|23|), all terms 
carrying suitable indices for the field they belong to, and where N\2 denotes an 
additional normalization factor. The latter is necessary, since we included an 
essential exponential describing the two-point correlations between the system 
and the environment fields. We obtain: 

N 12 (t) = (det{l - G!(t)C? 12 (i)G 2 (t)G 12 (*)}) 1/4 , (27) 

similarly as in eq. ( [25|) . 

Employing eq. (|27|), the calculation of the two-field effective action T is 
straightforward, even if tedious. We present here only the result for the sim- 
plest bilinear interaction, V($>\, $ 2 ) = — /^ 2 $i < I ) 2 , while the case of a general 
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quartic interaction is studied in RefO(a): 



r[* i2 ] 




n $ - in 2 - - 

LL 3 ^] 2 LL J 2 




+ - tr [ tuGiGiGu } - Tiiihx [G x G 2 G l2 



(28) 



with vf ] = d n v j (<f> j )/d<f>J l , (/} = Jd d xf{x)/V d , and V d = Jd d x. We 
also use the abbreviations G\i = G\2[A — B] , Gj = [A — S] _1 Gj , and V 2 
acts on either one of the two formal spatial coordinates of Gj ; furthermore, 
Si = S1-IG2G12S12 , and £ 2 follows by 1 ^ 2 ; finally, A—B = 1-G 1 G 2 G? 1 2 2 
and ^4 + B = 1 + G1G2S12 ■ We recall that 'products' of two-point functions 
involve integrations over intermediate coordinates, which we suppressed as 
before. The two-point functions are assumed to be translation invariant (bulk 
matter) . 

Even for the simple bilinear interaction the effective action is quite com- 
plicated due to the full Hartree-Fock approximation for the system and envi- 
ronment fields. Note the 0(K) and 0(h 2 ) quantum corrections to the classical 
action, which appears in the first line of eq. (|2^). We observe that the dress- 
ing of the two-point functions G\fii by a geometric series of terms involving 
themselves and G12 disappears, as soon as the latter correlation function van- 
ishes. The field equations resulting from the variations of the effective action 
are equally involved. In Section 4.3 we present a simplified version of them 
and some intriguing numerical results concerning the zero-dimensional limit, 
i.e. the quantum mechanical Hartree-Fock approximation. 

4-2 The Entropy Functional 

Knowing formally the wave functional of the system ($1) coupled to the envi- 
ronment ($2), we presently evaluate the von Neumann entropy of the system 
in terms of the two-point functions of the previous section and according to 
the outline in Section 2 . 

First of all, the system density functional is obtained from the total density 
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functional (matrix) by tracing over the environment d.o.f., cf. eq. (^): 

p s [$i,$i;t]= y2?* 2 (|p(t)|$i,*2> = y *t,[#i,# a ;*]*ia[*i,*a;*] , 

(29) 

which is a Gaussian integral again. We obtain: 

Ps[*i,*i;t]=*2 fl [«i;t]*G 1 [*i;*] ex P {r 1 *[$' 1 ;t]G 2 (t)y 1 [$ 1 ;t]} , (30) 
with: 

Y 1 [fr,t] = ±[*-* 1 ][Gi 3 (t)-iH 13 (t)] , (31) 



and where the effective Gaussian vp Gi here is defined as before, cf. eqs. (23) 
and (p^), however, with the replacements: 

JVi(t) — ► JVx(i) = iVi(t)JV 12 (t) , 

Ei(t) — Ei(t) , (32) 



with A and Ei as introduced after eq. (|28|). Note that the result of eq. fl3l| ) 
has the typical form of a modified pure state density matrix times an expo- 
nential influence functional; both modifications vanish in the limit of vanishing 
correlations between system and environment. 

In order to evaluate the entropy, i.e. -Trp^lnpg , we need to diagonalize 
the density matrix, such that its eigenvalues are accessible. This calculation 
was performed by functional means in Ref.H(b) . The result can be represented 
in the form: 

^) = -/^{Ml-n) + I ^lnn} , (33) 
where Xk denotes the d-dimensional Fourier transform of X(x), and we find: 

y * - wT-Btfn (34 > 

« K [G 1 G 2 (G 1 2 2 +E 1 2 2 )] fe , (35) 

in terms of A,B as before, and where the constant n is 1/4 (1/2) in the 
small (large) entropy limit. We employed here that the products of translation 
invariant two-point functions involving integrations over intermediate coordi- 
nates become convolutions, which factorize after Fourier transformation. The 
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result in eq. ( |33| ) represents a sum-xif oscillator like terms, which could have 
been expected, cf. the first of Refs.E.3. 
Several remarks are in order here: 

• The results obtained are formally independent of the dynamics, which 
only enters through the actual time dependence of the two-point func- 
tions (the 'mean fields' do not contribute); however, their gener- 
ality is restricted by the underlying ansatz of a Gaussian wave functional 



(TDHF), eq.(27) 



For vanishing correlations between system and environment (independent 
subsystems), G12 = £12 = , we obtain Ss{t) = . 

For small widths of the Gaussians, Gi or G2 small, the entropy is small. 
In this case the system or the environment follows a quasi-classical tra- 
jectory in 'coordinate' (field) space; the widths cannot be squeezed to 
zero because of the uncertainty principle, which is incorporated TDHF. 

Employing the entropy, eq. (|34|), we define a dynamical time scale: 

T -l d 1n o , /W /Wfc fY k /(l-Y k ) 

which determines the decoherence time td and the equilibration time r equ in 
the limits of small and large^entropies, respectively, as indicated. 

The calculation of Ref.tl(b) yields as a side product the most probable 
eigenstate of the system density matrix, i.e. the field pointer state with the 
largest probability and the smallest (field) kinetic energy. It is a coherent state, 
i.e. Gaussian wave functional centered around the classical field configuration 
<&i(x,t) with momentum Hi and effective real width: 

Geff = 1 - GiGatGa - 4)- (G!G 2 G 12 E 12 )2 ' (37) 

to be compared with the bare width 4Gi , in the absence of correlations. The 
special role of coherent states as pointer states of the electromagnetic field has 
also be found more recently in RefEj, in a model of a dielectric medium. 



4-3 Equations of Motion and Semiquantum Chaos 

In order to illustrate the dynamical content of the effective action derived in 
Section 4.1 , we consider simplified versions of the resulting equations of motion 
in this section, in particular also the case without environment. The full set 



1G 



of equations for arbitrary quartic interactions within and between system and 
environment can be found in RefO(a). 

Presently, we make the additional assumptions: i) We study only the in- 
frared limit, i.e. quantum mechanics of spatially homogeneous coupled fields 
$1, $2 ■ h) We neglect mean fields in the environment, $2, II2 ~ 0, and con- 
sider the case of small correlations G12, £12 between system and environment. 

Varying the effective action, eq. (|8|) , and applying the above assumptions, 
we find the set of coupled nonlinear (first order) equations of motion : 



ST 



=^ Mi=Hi , (38) 



§-=0 =► ft* = - \v^ Gl , (39) 



ST 



= => d t Gj = AEjGj , (40) 



ST 1 
-—=0 => d t Gi2 = -2(£ 1 + £ 2 )Gi2--(G 1 - 1 +G 2 - 1 )£i2 , (41) 

Furthermore (j' ^ j, j = 1,2): 
^ = =► ftE, = -2S/ + Igt^ _ ^f) _ Gj - ^G f G 12 , (42) 

<sr 1 

— - = => ftE 12 = -2(Ei + S 2 )S 12 + -(Gr 1 + Ga x )Ei2 + 2 M 2 , (43) 

0Lr\2 & 

where all quantities are simply functions of time; we combined the resulting 
equations in such a way that time derivatives on the r.h.s. were eliminated. 
Note that the terms oc /i 2 are due to the system-environment coupling. 



Several features of eqs. (38)-(|43|) appear to be of a rather general nature. 
First of all, the 'momenta' fti, Ej can be easily eliminated. The resulting 
second order equations are of the (inhomogeneous) anharmonic oscillator type. 
They potentially become unstable due to a dynamical 'negative mass squared' 
term; in this case the effective potential corresponding to T has a complicated 
structure allowing for tunneling processes E2l. 

Furthermore, the coupled linear eqs. (§d|), ( [43| ) can be integrated formally. 
We consider here only the solution of one of them in terms of the other: 

Gl2 (t) = G 12 (0)e- 2 /o dt"(^+^)_i jt dt , Eia (Gr i +G ,-i )c -2/;dt»( Sl+Sa)j 

(44) 
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which shows that the evolution of the system-environment correlations, and 
therewith of the whole dynamics, involves characteristic non-Markovian fea- 
tures. Furthermore, one may expect a large number of nonzero components in 
the Fourier spectra of the two-point functions of the system or environment 
alone. The functions G12, £12 inherit these and, therefore, the non-Markovian 
behavior leads to a quasi stochastic influence of the environment. 

Finally, expanding G12, £12 for sufficiently short times, we explicitly solve 
eqs. (fll|), (f43|). In the limit of small but finite correlations we obtain: 

G 12 (t) « G 1 2-2[(£ 1 +E 2 )G 12 + i(Gr 1 +G 2 - 1 )E 12 ]i 

-i M 2 [Gr 1 + G 2 - 1 + 0(G 12 ,I] 12 )]t 2 , (45) 

E 12 (t) « £ 12 + 2[ A1 2 -(£ 1 + £ 2 )£ 12 + i(G- 1 +G 2 - 1 )G 12 ]i 

-2/i 2 [E 1 + S 2 + 0(G 12 ,I] 12 )]t 2 , (46) 

where all two-point functions on the r.h.s. here assume their initial value 
(t = 0). Note how the analytical behavior changes as Gi 2 (0), Ei 2 (0) — > . 

This is reflected in the decoherence time tq of our simple model, which we 
calculate employing eqs. (|34|)-(|36|) of the previous section and assuming small 
but finite Gi 2 , Ei 2 : 

r -i « ^-ln[GiG 2 (G 12 + E»)]m. « 4(Er(0) + £ 2 (0)) + 4 /'-^ (,)1 



dt ^^ -v ^ u , ~ la/J ™ ■ ~,v-yy . Gl2 (0) 2 + Sl2 (Q) 2 ■ 

(47) 

In the absence of the initial correlations Gi 2 (0), Ei 2 (0) the decoherence time 
vanishes instead, td ~ t/2 (t — ► 0). This demonstrates the important role of 
correlations, which has been investigated in more detail in models of quantum 
Brownian motion (cf. Section 3), see e.g. Refs.tHlO. For an environment 
without self- interaction at high temperature one finds td oc TEJ. 

For the remaining part of this section we neglect the environment, in order 
to point out an important effect of the TDHF approximation employed in 
the derivation of the equations of motion (|38"|)-([43"|). Considering the effective 
action, eq. (p8|), it has become obvious that this approximation includes only 
the 0(K) and 0(h 2 ) corrections to the classical action. 

Whereas the classical equations of motion of the system alone (Id anhar- 
monic oscillator = 2 d.o.f.) £an have only regular solutions according to the 
Poincare-Bendixson theoremtj, the potential for scmiquantum chaos is obvious 
in the case of the four coupled nonlinear first order equations, cf. eqs. (^8|)-(^0|) 
and (|42]) , representing the systenoJn TDHF. A detailed study of this subset of 
equations has been made in Ref.c2l. 
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Here we represent some illustrating numerical results obtained by inte- 
grating the equations of motion forward in time. The initial configuration is 
a Gaussian wave packet with conserved energy E (all quantities are suitably 
rescaled and presented in dimensionless form). For the chosen model parame- 
ters the minimum of the energy is at E = —24.3 , the top of the hill separating 
the two wells of the classical double-well oscillator potential is at E = 0.0 . 

In Fig. 4 Poincare sections for various energies are shown (0 = $(t), 7r = 
ft(t)). We observe the characteristic break-up of KAM tori as the energy_is 
increased and the corresponding stochastic filling of the 4>, tt phase space 
This clearly shows a transition between regular and chaotic motion, which can 
be associated with tunneling paths, despite the fact that the present classical 
model is regular E3. 

Furthermore, in Fig. 5 the behavior the Fourier transform w.r.t. time of 
a diagonal element of the density matrix is shown. At very low energies (not 
shown) one finds a simple line spectrum involving the fundamental frequency 
ujo (i.e. harmonic approximation of the potential) or a few multiples thereof 
(regular motion). At the energies shown the full trajectories are chaotic and 
produce a broadband noise spectrum. At very high energies the motion be- 
comes regular again with the wave packet experiencing essentially only the 
anharmonic quartic part of the potential. 



E=-12.49 E=-6.61 




4> 4> 

Figure 3: Poincare sections for the system at various energies (after Ref. [20]) showing the 
transition to semiquantum chaos, cf. main text. 
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Figure 4: Fourier transform of a diagonal element of the system density matrix at various 
energies (after Ref. [20]) showing semiquantum chaos as broadband noise in the expected 
line spectrum (nonlinear resonance at E=-3. 68601), cf. main text. 



These results raise the question of the reliability of the underlying TDHF 
approximation or analogous semiclassical perturbative expansions around the 
bare (h = 0) classical limit. Generally, the full quantum evolution according 
to the Schrodinger equation incorporating an environment implies a linear but 
non-Markovian master equation for the system density matrix, which might 
rule out the possibility of 'quantum chaos'. 

Furthermore, it has recently been shown that any variational approxi- 
mation for a closed system is equivalent to. an effective classical Hamiltonian 
dynamics for the variational parameters 123. There, as we have seen in the 
TDHF case, eqs. (|38|)-([43|), the potential for semiquantum chaos arises, if the 
underlying classical Lagrangian involves higher than quadratic terms. In some 
examples then, the onset of semiquantum chaos has indeed been related to the 
breakdown of semiclassical approximationsEa, even if they do not require weak 
coupling. 

However, the important question remainsE§, whether there are situations 
involving (infinitely) many d.o.f., where such approximations could reliably 
represent the highly complex, even if overall quasi-periodic quantum evolution 
of the coupled system+environment (fields). Such cases have recently been 
found in solid state devices. Nevertheless, the real time evolution of fields 
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which are relevant for high-energy experiments still needs more study, in order 
to understand decoherence processes, entropy production, and thcrmalization 
quantitatively. 

5 Conclusions 

We have seen in Sections 1,2 that entropy production in a quantum system 
can be understood dynamically on the basis of von Neumann's definition of the 
entropy in terms of its density operator, if and only if the system is coupled 
to an environment of unobserved or even unobservable degrees of freedom. It 
induces the necessary quantum decoherence EJ. 

The example of a nonrelativistic quark in Brownian motion in a simple 
model of a gluonic environment has been presented (Section 3). It has been 
pointed out that an environment with a spectral density distribution which 
is enhanced in the infrared, as compared to the frequently studied 'Ohmic', 
phonon, or photon environments E3, may lead to interesting squeezing, stabi- 
lizing, or enhanced spreading of the wave packet effects, which affect the rate 
of entropy production. 

In order to extent the decoherence approach to field theory, a system and 
environment consisting of two distinct scalar fields have been studied in the 
time dependent Hartree-Fock approximation in Section 4u. We derived the 
effective action including up to 0(h 2 ) corrections, the equations of motion, 
and the entropy functional and discussed their general properties. 

Particular attention has been paid to the fact that the approximate semi- 
classical dynamics may be intrinsically unstable and leading to deterministic 
semiquantum chaos. It might signal a breakdown of any such approxima- 
tion scheme for genuinely nonlinear interacting field theories. Whether under 
these circumstances the initial decoherence, entropy generation, and possibly 
thermalization in high-energy hadronic or nuclear collisions can be reliably 
calculated, is an interesting question. The decoherence effects due to initial 
bremsstrahlung mentioned in Section 2 are presently studied. 
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